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Abstract 

The multiple passages of light through an absorption inhomogeneity of finite 
size deep within a turbid medium is analyzed for optical imaging using the 
"self-energy" diagram. The nonlinear correction becomes more important for 
an inhomogeneity of a larger size and with greater contrast in absorption with 
respect to the host background. The nonlinear correction factor agrees well 
with that from Monte Carlo simulations for CW light. The correction is about 
50% — 75% in near infrared for an absorption inhomogeneity with the typical 
optical properties found in tissues and of size of five times the transport mean 
free path. 

(c) 2008 Optical Society of America 
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The main objective of optical imaging of turbid media is to locate and identify the 
embedded inhomogeneities by essentially inverting the difference in photon transmittance 
in time or frequency domains due to the presence of these inhomogeneities . ^i^i^i"^ The key 
quantity involved is the Jacobian which quantifies the influence on the detected signal due 
to the change of the optical parameters of the medium. The linear perturbation approach is 
only suitable to calculate the Jacobian for a small and weak absorption inhomogeneity, and 
not valid when the absorption strength is large.— This failure can be attributed to the multiple 
passages through the abnormal site by the photon. The most important correction is the 
"self-energy" correction,— which takes into account the repeated visits made by a photon 
through the site up to an infinite number of times. The presence of other inhomogeneity 
"islands" can be ignored because the photon propagator decreases rapidly with the distance 
between two separate sites. 

In this Letter, the nonlinear correction for an absorption inhomogeneity of a large strength 
due to repeated visits by the photon is modeled by a nonlinear correction factor (NCF) to 
the linear perturbation approach. NCF as a function of the size and the strength of the 
inhomogeneity is estimated using the "self-energy" diagram. NCF is obtained based on the 
cumulant approximation to radiative transfer, and verified by Monte Carlo simulations for 
CW light. The magnitude of NCF is 0.5 — 1 for an absorptive inhomogeneity of size up to 5/f 
{It is the mean transport free path of light), and of the typical optical properties of human 
tissues {fiJt/c ~ 0.01 — 0.05 where Ha is the absorption coefficient and c is the speed of light 
in the medium). 

Consider an absorption site centered at r and far away from both the source and the 
detector, the change of the detected light. A/, at the detector from a modulated point 
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source at including the multiple passages through the site, is given by 



AI = -G(r,,a;|f)l^5^,(f)£ [-iVseif(cu;i?)l^<5/i,(r)]"G(f,a;|r,) 



(1) 



n=0 



1 + Nseu{^;R)V6fia{r) 



where (5/ia is the excess absorption of the absorption site of size R and volume V, uj is the 
modulation frequency of light, G is the propagator of photon migration in the background 
medium, and 



is the self-propagator which describes the probability that a photon revisits the volume V. 
Here G'(r2, c^lri) gives the probability density that a photon leaves the volume at ri and 
reenters it at r2. The scattering property of the site is the same as that of the background. 
In Eq. (P) G{rj^, uj\f) and G{f,uj\rs) are well modeled by the center-moved diffusion model as 
long as the separations jr^ — f | , jr^ — r| ^ /i.- However, the diffusion Green's function can 
not be used in Eq. (j2I) to evaluate iVseif(ci;;-R) because the diffusion approximation breaks 
down when ri is in the proximity of r2. 

Comparing Eq. (0) to the standard linear perturbation approach, the nonlinear multiple 
passage effect of an absorption site is represented by a nonlinear correction factor (NCF) 



This factor serves as a universal measure of the nonlinear multiple passage effect as long 
as the absorption site is far away from both the source and the detector and its size is 
much smaller than its distance to both the source and the detector. This correction is more 
significant when NCF is further away from unity. 

The photon propagator A^(r2, t|ri, s), the probability that a photon propagates from 
position ri with propagation direction s to position r2 in time t, for any separation between 




(2) 



NCF= l + N,,u{u;;R)V6fia{f) 



(3) 
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ri and r2, was recently derived^ in a form of the cumulant approximation to radiative 
transfer. 

In the case of interest where the absorption site is deep inside the medium, the photon 
distribution is isotropic. The photon propagator is simphfied to N^s^r, t) = Nca{\'Y2 — ri| , t), 
obtained by averaging A^(r2, t|ri, s) over the propagation direction s of hght over the 47r sohd 
angle. In the frequency domain, this effective propagator is approximately given by 

, , 4^exp(-i/€V) + ^lfe^sinh(/cr) r </t 

where D = ltc/3 and k = [3(/ia — i^^)/kc]^^'^ whose sign is chosen with a nonnegative real 
part. The two terms in N^s when r < U represent ballistic and diffusion contributions 
respectively. The ballistic term does not depend on scattering because the photon distribu- 
tion involved is already isotropic. Only diffusion contributes to N^^ when r > If. The self 
propagator for an absorption sphere deep inside the medium is given by: 

N,,ii{u;R) = ^ J^J^N,s{\r2-ri\,u)dh2dhi 

= TT / Ne{{{r,uj)-fo{r)^7rr^dr (5) 
V Jo 



where 7o(r) = + 3^ (■^) is the characteristic function for a uniform sphere^. An absorp- 
tion site of an arbitrary shape can be treated in the same fashion. The exact self propagator 
must be computed by a numerical quadrature. A good approximation of iVgeif (ci^; R) is 

NM^-R) = -^\ > (6) 

i (1^^ + I - + U-'') - + 0{.^) ^ > 1/2 

using Eq. (@)) where ^ = R/lt when \k,\ <^ 1. This approximation (jUI) compares favorably 
to the exact iVseif(co' = 0]R) obtained by numerical quadrature for continuous wave light 
propagating inside a nonabsorbing medium. The exact and approximate versions of the 



dimensionless self-propagator NseuVlt^c are plotted as solid and dashed lines, respectively, in 
Fig. ((T^). The dimensionless self-propagator iVgeifV^/t^^c depends solely on two dimensionless 
quantities ult of the background and R/lt of the absorbing sphere. 

It is worthwhile to point out that the self-propagator in time Nseif{t]R), the inverse 
Fourier transform of includes the contribution from the ballistic motion of the photon 
when the photon passes through the site. This ballistic contribution manifests itself as 
the linear decay of Ngesit; R)V in the form of 7o(ct) near the origin of time, followed by a 
transition to diffusion [Fig. (db)]. 

The nonlinear correction factor is obtained by plugging (jSj) or (jH)) into In particular, 
we have 

NCF^ L i6.n^ s^J ? - / 

[ [i + (r^ + - + U-'T ^ > V2 

where q = VS^aij)/^^ is the dimensionless strength of the absorber when |k| i? -C 1. The 
effectiveness of an absorber of a fixed strength depends on its size. For an absorber of a fixed 
g > 0, the effectiveness of absorbing light is diminished (NCF decreases) when its size is 
reduced. This can be understood from the fact that the photon spends less time per volume 
inside the absorber of a smaller dimension because of the ballistic motion of the photon 
after each scattering event. The photon leaves a small site {R < k) in an almost straight 
line. The diffusion behavior for an individual photon is only observed after a large number 
of scattering and on a scale larger than It. 

Fig. Q shows plots of NCF versus absorber size for typical absorbers of excess absorption 
^fJ'Jt/c equal 0.01 and 0.05, respectively. The nonlinear correction factor generally decreases 
with the size of the absorber whose excess absorption is fixed. With the increase of the 
background absorption and the modulation frequency, the nonlinear correction becomes less 
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accentuated. The phase delay is larger for higher modulation frequencies and less background 
absorption. 

Monte Carlo simulationsi^ are performed for CW light propagating in a uniform nonab- 
sorbing and isotropic scattering slab. The thickness of the slab is L = 80/^. A spherical 
absorber of radius R is located at the center (0, 0, L/2) of the slab. The excess absorption of 
the absorber is 6fiJt/c = 0.01. The absorber has the same scattering property as the back- 
ground. The details of the MC computation has been provided in a previous publicationii. 
The correlated sampling method is used in each simulation to reduce variance.— A single 
simulation is used to compute the change in light transmittance due to the presence of the 
absorption site and the corresponding nonlinear correction factor. Fig. (jSt^) shows the non- 
linear correction factors obtained from numerical quadrature, the approximate form Eq. 
and Monte Carlo simulations, respectively. The agreement between our theoretical NCF 
and that from Monte Carlo simulations is excellent. The slight difference between them at 
large radii is accounted for by the fact that the sphere can no longer be regarded as small 
compared to the dimensions of the slab. The probability of a photon revisiting a large sphere 
is overestimated by Eq. (j^ for the sphere located at the center of the slab.— 

Fig. (jnb) shows the percentage change of the CW transmittance estimated from the 
experimental data given in Fig. (9) of Ref . (j5|) . The relevant parameters of the experiment are 
summarized in the inset. The theoretical predictions from the linear perturbation approach 
with and without the nonlinear correction are also shown in Fig. (j^b) assuming a coUimated 
point source and a point detector in a cofocal setup. The agreement between our theoretical 
prediction with nonlinear correction and the experimental result is good. 

The typical value of the absorption coefficient of human tissues in the near infrared 
indicates that fiak/c ~ 0.01 — O.OS .^^i^^ This fact should put our results on NCF in this 



range [Figs. (j21) and Q] into perspective. Our analysis reveals that the nonlinear correction 
for absorption inhomogeneities due to multiple passages of light in optical imaging of human 
tissues becomes important when the size of the inhomogeneity is much larger than the mean 
transport free path of light. This factor is smaller (further away from unity) for higher excess 
absorption and weaker background absorption. The value of NCF decreases from ~ 0.75 
to ~ 0.5 for an absorption site of radius 5lt with excess absorption 6fiah/c increasing from 
0.01 to 0.05. This corresponds to 75% and 50% underestimations of the excess absorption 
coefficient in these cases if the standard linear perturbation approach is applied naively. 

In conclusion, the nonlinear correction is analyzed for the effect of the multiple passages 
through an absorption inhomogeneity of finite size deep inside a turbid medium in optical 
imaging using the "self-energy" diagram. The nonlinear correction factor is verified by 
Monte Carlo simulations for CW light and supported by experimental results. The nonlinear 
correction becomes more important for an inhomogeneity of a larger size and with greater 
contrast in absorption with respect to the background. The standard linear perturbation 
approach in optical imaging should be augmented to include this nonlinear correction. 

This work is supported in part by NASA and Army. M. X. thanks the support by the 
Department of Army (Grants DAMD17-02- 1-0516). The U. S. Army Medical Research 
Acquisition Activity, 820 Chandler Street, Fort Detrick MD 21702-5014 is the awarding and 
administering acquisition office. The authors are indebted to the anonymous referees who 
help in improving this paper. 
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Fig. 1. (a) The self propagator Nseif{uj; R)Vl^^c and its approximation form when k — 0. 
(b) The self-propagator for spheres of various radii in the time domain inside a nonabsorbing 
medium. 
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Fig. 2. The nonlinear correction factor (magnitude and phase angle) vs the size of absorbers 
whose excess absorption (5//aZt/c equals 0.01 and 0.05, respectively. Note k^I^ — 3(/Ja— ia;)Zt/c 
for the background medium. 
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Fig. 3. (a) The theoretical nonhnear correction factors from numerical quadrature (Exact), 
the approximate form Eq. ((7j) (Approx), and Monte Carlo simulations (MC). Results from 
four independent MC simulations are shown for each radius. The standard linear perturba- 
tion approach corresponds to a horizontal line NCF = 1 (not shown in the figure), (b) The 
percentage change of the CW transmittance from the experimental data given in Fig. (9) of 
Ref. (p) compared to the theoretical predictions made by the standard linear perturbation 
approach (StdPert) and the ones including NCF (Exact and Approx). 
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